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^ . ABSTRACT 

^ ; 

Context. The metal-rich bulge globular cluster NGC 6441 shows a well developed blue horizontal branch (Rich et al.), together with a strong 
slope upward from the red clump to the blue of the RR Lyrae region. Both features, the former corresponding to the well-known second 
parameter problem, are not explained by conventional evolutionary models. 

Aims. Helium self-enrichment is proposed as a possible solution to both questions, a mechanism already invoked for the interpretation of the 
^ peculiarities in NGC 2808 and M 1 3 . 

Methods. We make use of horizontal branch simulations, covering a wide range in main sequence helium abundance, to investigate whether 
the main features of NGC 6441 horizontal branch population, including the RR Lyrae variables period, can be reproduced. 
^ I Results. To describe the horizontal branch of NGC 6441, the helium content Y in the red clump must reach at least 0.35; values up to Y~0.37 
' are necessary to populate the RR Lyr region, reproducing also the observed mean long periods; depending on the dispersion in mass loss 
t| assumed in the simulations, values up to Y~0. 38-0.40 are necessary to populate the blue HB. The total self-enriched population amounts to 
~60% of the whole stellar content. 

, Conclusions. Self-enrichment and multiple star formation episodes in the early evolution of globular clusters appear more and more able 
H* ■ to account for many of the chemical and population peculiarities observed in these systems. The very large helium abundances (Y t 0.35) 
. 1^ required for ~ 14% of the horizontal branch population pose some problem on the enrichment mechanisms. 

X 
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1. Introduction 

In the last few years, some of the (by now) well known abun- 
dance anomalies in globular cluster (GC) star composition have 
been observed in unevolved cluster members (f.e., the anticor- 
relation between Na and O abundances; Gratton et al. 2001, 
Ramirez & Cohen 2003, Cohen & Melendez 2005, Carretta 
et al. 2006). These observations have focused attention on the 
early stages of GC formation and evolution, since it appears 
necessary that the abundance deviations from that of the pris- 
tine material have taken place before the formation of the stars 
in which they are observed. 

A model which can explain all the chemical anomalies is 
not yet available, but some aspects of the problem begin to be 
understood, while others are still obscure. A clear step forward 
was made when it was recognized the link between the chem- 
ical anomalies and the horizontal branch (HB) morphology of 
the clusters: if the chemical anomalies, and in particular the 
Na-O anticorrelation, is interpreted in terms of a second gen- 
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eration of stars born from the ejecta of the massive asymptotic 
giant branch (AGB) stars of the first stellar generation, there 
must also be a helium discontinuity between the two samples 
of stars. In fact the AGB ejecta from which the second gener- 
ation is born are all helium enriched 1 to'Antona et al.U2002h 
and consequently the second star generation evolving along the 
HB have bluer locations with respect to the first generation 
HB stars. This model was able to produce a detailed expla- 
nation of th e anomalous mor phology of the HB of the cluster 
NGC 2808 ( lD'Antona& Cal oi. 2004); besides, observations 
have confirmed the existence of a Na-O anticorrelation in the 
cluster (Carretta et al,. 2006) . 

An intriguing problem concerning at the moment only 
NGC 2808 and the "cluster" co Cen is the presence of a well 



' The masses of the AGBs involved must be i 4M0, otherwhise the 
second generation would show evolutionary signatures of a prolonged 
phase of third dredge up, namely the sum of CNO abundances would 
not be costant, and there would be s- process elements enhancement. 
The helium content of the ejecta of a 4 Mg AGB star is Y~0.27 for 
Z=0.001. 
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defined second main sequence (MS), bluer than the main one. 
In terms of standard MS models, this feature is a helium rich se- 
quence in which the helium mass fraction must be Y~0.4 (see 
for 0) Ce n Norris 2004; Piotto et al. 2005 and for NGC 2808 
Id 'Antona et alJ2005l) . AGB models do not predict such a large 
helium abundance in the ejecta, the maximum computed value 
being Y^O.36 in Lattanzio et al. ( 2004), obtained from the evo- 
lution of a 6M0 with Z=0.004. Ventura et al. (2002) quote val- 
ues up to Y^0.31 and Ventura & D' Antona (2005b) estimate 
Y=0.32 for their most massive models (6.5 Mq) for metallic- 
ity Z=0.001. Although this value can be considered conserva- 
tive, as these models do not include any kind of overshooting, 
a preliminary exploration of the second dredge up does not in- 
crease Y beyond ~ 0.35 (Ventura 2006, private communica- 
tion). The episodes of third dredge up may help to increase the 
helium content in some models (Ventura & D' Antona, 2005a), 
but too many episodes in the end do not preserve the constancy 
of C-hN-i-O. Inclusion of the results from these computation in a 
code for computing the chemical evolution of helium confirms 
the conclusion that present A GB models cannot ex plain the ex- 
treme helium enhancements jKarakas et al.Ll2006 l). The role of 
super-AGBs, the stars which ignite Carbon in semidegenerate 
conditions and may evolve through a AGB phase, ending up 
their life as O-Ne-Mg white dwarfs, has not yet been explored 
for the metallicities of GCs. Siess (2006) finds Y ~ 0.38 for 
these stars after the second dredge up, but this result may be af- 
fected by the initial value of the helium content, as his models 
refer to s olar metalhcity. 

From lP' Antona et alJ (|2005), the extreme helium enhance- 
ment, derived both from the blue MS and from the extreme HB 
blue tails EBT2 and EBT3 (as defined in Bedin et al. 2000), 
concerns only ~15% of the stars in NGC 2808, and a similar 
fraction of stars is involved in the bluer MS in oj Cen. 

We remind that the extreme value of the helium content, 
obtained from the MS fit, is rather uncertain. On the basis of 
our present understanding, extreme helium enhancements (Y > 
0.35) may derive either from the high mass tail of AGBs (for 
which at present we do not yet have models, see above) or from 
the wind ejecta of the massive stars of the first stellar genera- 
tion (Piotto et al., 2005; Prantzos & Charbonnel, 2006; SmitlJ 
I2OO6 ). This latter suggestion would complicate even more the 
picture of star formation in the "simple" stellar population such 
as GC stars were considered until a few years ago. It is then 
important to ascertain whether there are any other GC in which 
there is a clear indication for an extreme helium rich (Y>0.35) 
component. In the absence of new information coming from 
other main sequences, we considered modelling the peculiar 
HB of the cluster NGC 6441, whose morphology and RR Lyr 
population (Rich et al. 1997, Layden et al 1999, Pritzl et al. 
2000) are both a clear indication of the presence of helium en- 
richment, as alr eady suggested bv^Sweig art & Catelan L1998) : 
ID' Antona et"al] (|2002); D' Antona & Caloil ( |2004 . The results 
obtained from this analysis can be extended to the cluster NGC 
6388, very similar to NGC 6441 in metallicity and peculiar HB 
morphology and RR Lyr periods. 

We shall proceed as follows: §2 presents the observational 
data for NGC 6441 and the related problems; the theoretical 
background is given in §3, while §4 describes the comparison 
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Fig. 1. CM diagram of NGC 6441 from the data by Piotto et al. 
2002. The crosses (in red) indicate our selection of HB mem- 
bers; also shown is the histogram of their distribution with re- 
spect to B-V 
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Fig. 2. An enlargement of the HB red clump region: crosses 
(red) indicate stars with B-V > 0.5 mag, open circles (blue) 
stars with B-V < 0.5 mag. The histogram of red HB stars with 
respect to V mag is shown. 



with the observations with a possible solution to the questions 
in §2; some aspects of the given solution are examined in §5, 
and §6 gives the final discussion. 

2. The horizontal branch in NGC 6441 

NGC 6441 belongs to the metal rich group of GCs (Z i 
0.004). All these clusters have a HB confined to the red side 
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of the RR Lyr region, with the exception of NGC 6441 itself 
and of its "twin cluster" NGC 6388 (Rich et al. 1997). 

The CM diagram of NGC 6441 is shown in Fig. [T] (data 
from Piotto et al. 2002, dereddenedby E(B-V)=0.44). The red 
clump has a strong slope upward, which continues in the RR 
Lyrae and blue HB regions, with a difference of at least 0.6 
mag between the top of the blue HB and the bottom of the red 
clump. The features to be explained are: a) the presence of HB 
stars hotter than the red clump ones; b) the strong tilt in lumi- 
nosity, much larger than expected on the basis of evolutionary 
models (f.e., Raimondo et al. 2002); c) the very long average 
period of the RR Lyr variables, unexpected on the basis of the 
metal content (Z > 0.006; Layden et al. 1999, Pritzl et al. 2001, 
2003); d) the presence of an extended blue tail down to a mag- 
nitude V ~ 19.5, corresponding to My ~ 3.5 mag. 

In Fig. n our choice of the HB members is indicated 
(crosses). The numbers turn out to be: a total of 1452 mem- 
bers, out of which 1278 form a red clump ((B-V)o > 0.5). 
For our simulations we assumed therefore a total of 1452 stars, 
to be distributed according to observations. Also in Fig.^the 
histogram of the number of HB stars with respect to B-V is 
shown; in Fig.|2lthe histogram with respect to V mag is shown, 
considering only stars with (B-V)() > 0.5 (crosses)^ 

Th e number of RR Lyr is 45 in the sample bv lPiotto et al.l 
ll2002h (Piotto, private communication). We assume a mean 
fundamentalized period (see below ) of 0.67 d, based on 63 
RRab and RRc bv lPritzl et all (l2003l) . 

3. The evolutionary and synthetic l-IB models 

The computation of HB models and of the synthetic HBs fol- 
lows the prescription in D'Antona & Caloi (2004). We will 
model the HB by assuming cluster members to be a mix- 
ture of stars with different helium abundances, as we did 
for NGC 2808 and Ml 3. Previous simulati ons for the HB 
of NG C 6441 have been performed by Sweigart & Catela*!] 
ill 9981) . who also considered models with high Y, but not a mix- 
ture of helium contents. 

The scheme of our modelling will be the following. We 
consider a first group of stars with a fixed helium abundance 
Y, close to the cosmological value Y= 0.25, which represents 
the first stellar generation. A second (or more) generation(s) 
are added, with a helium content variable from Y i 0.27 to 
Y i 0.40. The number vs. helium distribution will be chosen 
in order to comply with the relevant HB features mentioned in 
Sect. 2 

3.1. HB structures 

We computed evolutionary sequences from the main sequence 
to the red giant tip for a heavy element content of Z = 0.006 
and Y = 0.24, 0.26, 0.30 and 0.40. Having an estimate of the 
helium core masses at the flash, we obtained HB evolutionary 
sequences with the appropriate core masses and atmospheric 
helium abundances (after the first dredge-up). The theoretical 
tracks we re transformed into Joh nson magnitudes B and V by 
means of iBessell. Castelh. & Plez. (.1998) . The clump colours 
correspond nicely to the colours of the transformed tracks. 



Nevertheless, we do not emphasize this correspondence as a 
proof of a good choice of metallicity, as we could have shifted 
the observed colors within the uncertainties of the reddening 
(±0.03 according to Layden et al. ( 1999)). Such an uncertainty 
has no influence upon our comparisons, based quantitatively 
on luminosity diff'erences and periods and only qualitatively on 
the color extension of the HB, which greatly exceeds the red- 
dening error. On this basis, HB simulations were computed in 
V, B-V, assuming a gaussian error on magnitudes and colours. 
A standard deviation of 0.05 mag has been chosen on the basis 
of the HB data from Piotto et al. (2002). 

3.2. The main parameters in the simulations 

We assumed as the basic parameter describing the HB distribu- 
tion the ratio R of the number of blue ((B-V)o < 0.5) HB mem- 
bers to the total number of HB stars. This ratio (174/1452 with 
the numbers above) is R = 0.12 + 0.01. The other constraints 
to be satisfied by the simulations are the agreement with the 
histogram in Fig.|2 the extension in magnitude of the blue tail 
and the mean RR Lyr period. 

With regard to the RR Lyr variables, for every simulation 
we evaluated the mean period of the artificia l RR Lyrae vari- 
ables by means o f ivan Albada&Bakei< ll 19731) formulation; the 
hmits of the RR Lyr region were taken from logTeff = 3.79 to 
3.875. This is a standard interval for the more metal poor GCs, 
and may not correspond exactly to the present conditions of 
high Z and high luminosity. The properties of the RR Lyraes 
in NGC 6441 appear extreme when compared with those of the 
whole family of GCs (see, f.e.. Fig. 10 in Pritzl et al. 2001), 
but still within the limits of current interpretation. So we are 
confident that the order of magnitude of the estimated period is 
adequate. The observational mean period is built by averaging 
the observed periods, taking care to "fundamentalize" the pe- 
riods of the RRc variables by adding 0.12 8 to the logarithm of 
their periods jvan Alhada & Bakeilll973h. 

4. The simulations and the comparison with the 
data 

Two features, one observational and the other theoretical, affect 
strongly the comparison between observations and simulations: 

1 . the fact that the RR Lyrae variables and the blue component 
in NGC 6441 are only a tail of the red clump distribution, 
covering slighly more than 10% of the total HB population; 

2. the behaviour in the HR diagram of HB models rich in he- 
lium. 

The evolutionary path of the helium rich models may be con- 
fined in the red clump region, or may perform long loops, 
reaching very blue colours, for a difference in total mass 
of the order of 0.02 M^: a well known fact, see, f.e., 
ISweigart & GrossI fl976,) . This is shown also in Fig. |3 Both 
these features conspire to vary noticeably the RR Lyr number 
and mean period, even by only varying the random extraction. 
Consequently, we will try to fit first of all the observed fraction 
of stars at B-V<0.5 mag, and will impose the extension of the 
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Fig. 3. Evolutionary tracks for a few helium rich HB models: 
the long loops that they perform cover large portions of the 
allowed HB region. The tracks end at Yc=0.10. 




Fig. 4. Simulations show the effect on the evolutionary tracks 
of an increasing helium content: open triangles (red) indicate 
stars with fixed Y=0.25 (red histogram), filled triangles (blue) 
stars with increasing maximum helium content (Y^ax) (blue 
histogram), as labeled in the four panels. The black histogram 
refers to the whole HB. Z=0.006 in all models. 

blue tail down to My > 3 mag. In any case, we will not accept 
for the RR Lyrs a mean period smaller than 0.66 d. 

AA. Age, helium distribution and mass loss 

We assumed a reference age of 11 Gyr Coeval populations 
with the same metallicity but different helium content will 



have turnoff masses inversely proportional to Y. Therefore, the 
progeny of the helium rich, lower mass evolving red giants 
will be located, on the average, at larger temperatures (due to 
the smaller evolving masses) and larger luminosity (due to the 
higher efficiency of the hydrogen shell burning) than the struc- 
tures with lower helium. To illustrate this basic fact, we show 
in Fig.|4]the behaviour of two groups of stars, in one of which 
the maximum helium content is progressively increased. As ex- 
pected, the evolutionary paths, confined within the red clump 
for lower Y, as it happens for the cosmological helium popula- 
tion, reach farther and farther into the bluer HB regions when 
increasing the maximum Y. 

For an age of 1 1 Gyr, we chose as mass loss during the red 
giant evolution AM=0.165 Mq with cr=0.02. With this amount 
it was easy to obtain the desired behaviour of the clump lumi- 
nosity, when varying the helium content in the helium enriched 
population(s). A larger mass loss would cause an excessive loss 
of clump members (increasing the fraction of stars at B-V<0.5 
mag beyond the allowed values), a lower one would have the 
opposite effect. Of course the value of 0.165 is related to 
the chosen age: for 13 Gyr, a mass loss of 0.13 would have 
given similar results. 

The faintest blue HB stars in Fig. [Q- about 3.5 mag less 
luminous than the RR Lyrs - correspond to a HB mass < 0.5 
Mq. For Y - 0.40 and an age of 1 1 Gyr the evolving red giant 
has a mass of ~ 0.72 Mq, which would require a mass loss > 
0.22 Mq, larger mass losses would be necessary for lower Y, up 
to 0.45 Mo for Y=0.25. We have therefore considered models 
with Y up to 0.40, but even at this large helium content the very 
faint magnitudes are statistically difficult to obtain. So we took 
into account the possibility of an increase in the mass loss rate 
with the decrease of the red giant mass, according to Reimers 
(1977) mass loss relation (see also Lee et al. 1994). But this 
is not an automatic solution to the problem; the progression in 
mass loss used in Lee et al. (1994) would not allow any HB 
progeny to red giants of 1 1 Gyr with the assumed Z, helium 
content 0.40 and a basic mass loss of 0.165 M© (for the case 
Y=0.25), so we assume a milder dependence of the mass loss 
on the evolving red giant mass. Another constraint is given by 
the necessity of keeping the required number of stars within the 
red clump (and with the observed distribution in luminosity), 
since, with increasing mass loss, the evolutionary tracks tend 
to shift to the blue (see Fig.|3}- 

In Fig. |5] and Fig. |6l we give an example of a synthetic 
HR diagram and of its detailed behaviour in the clump re- 
gion. In addition to the parameters specified above — age 11 
Gyr, AM=0.165 Mq, cr=0.02 — , the extra mass loss included 
takes the form 0.10 x [Mrg(Y = 0.25) - Mrg(Y)]. Here we dis- 
cuss the main features of these diagrams. 

To reproduce the strong upward tilt in the red clump, it is 
necessary to consider objects with helium content up to Y ~ 
0.35. In Fig.[TOlthe typical relation of B-V with Y is shown. 
To get B-V < 0.5 mag, one needs stars with Y > 0.33, while 
the main body of clump stars (B-V > 0.5) have a Y content up 
to 0.37; to find (almost) exclusively blue stars we have to attain 
Y >0.37. 

In order to have, on the one side, blue HB stars at least ~ 
3 mag less luminous than the HB luminosity peak, and on the 
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other, the expected number of red clump members, it is nec- 
essary to assume a non uniform Y distribution for the second 
star generation. Otherwise, it is not possible to obtain a blue 
tail without remaining with too few red stars (B-V > 0.5). 

The difference in helium content between the first and the 
second generations has to be ~ 0.02 (a gap in Y from 0.25 to 
0.27), in order not to increase the peak in the clump distribution 
(Fig.0. At the same time, the gap cannot be larger than 0.03 
(Y from 0.25 to 0.28), otherwise the clump distribution would 
show (always) two well defined peaks (Fig.|8|i. 

The importance of the assumed error on V ,B-V is exem- 
plified in Fig.|5J where the same simulation of Fig.|6lis shown 
with an error on V ,B-V reduced from 0.05 to 0.03 mag: the 
clump appears clearly separated in two groups. This illustrates 
the relevance of improving the photometric precision in order 
to understand better the population stratification in GCs. 

4.2. Metallicity 

The adopted metallicity Z=0.006 falls on the lower side 
of the values ac ceptable for the cluster under exam (e.g. 
Iciemen^i ni et al.L 2005i). I n fact , it corresponds to the value 
in the ICarretta & Grattonl ( 1 19971) scale, if we do not consider 
a-enhancement. On the other hand, this metal rich cluster 
seems to be also a-enhanced (Carretta, private communica- 
tion). Should Z be larger, all the estimated helium values would 
be pushed towards even larger values. 



4.3. "The best solution" 

We consider a satisfactory result the one given by the simu- 
lation in Fig. Is] and Fig.|6l The helium distribution in shown 
in the lower panel in Fig. |6j an extra mass loss of 0.10 x 
[Mrg{Y = 0.25) - Mrg(Y)] has been assumed, beside the stan- 
dard AM=0.165 Mq mentioned before. The simulation has the 
correct value of the fraction of stars at B-V < 0.5 mag within 
the errors, a number of RR Lyr variables close to the estimated 
one (45 in the sample in Piotto et al., as mentioned before) with 
a mean period of 0.69 d close to the observed fundamentalized 
period of 0.67 d. As for the very blue stars, three stars are found 
at Mv > 3.0 mag in Fig. 13 

We stress again that variables and BHB stars are only a tail 
of the cluster star distribution, and that a small fluctuation in the 
clump number (± yfiNdump) and in the number of stars which 
leave the clump at a luminosity slightly lower than the average, 
changes noticeably the number of RR Lyraes and their mean 
period. 

4.4. Dispersion in mass loss and maximum helium 
content 

In the simulations we ado pt cr = 0.02 a widely used 
value in the relevant literature jLee. Demaraue. & Zina(.1994..) : 
ICatelan et al.l (1200 ih V With this assumption, it is necessary to 
include stars with Y up to 0.40. It is important to see whether 
such extreme values of Y can be avoided, since current models 
for possible self-enrichment in helium in globular clusters do 




Fig. 5. A synthetic CM diagram of the HB in NGC 6441, for 
an age of 1 IGyr, mass loss of O.165M0with cr - 0.02 and extra 
mass loss as in the text. Stars with fixed Y=0.25 (open trian- 
gles, red) plus stars with variable Y up to 0.40 (stars, blue). For 
the helium distribution, see text and the following figures. The 
standard deviation assumed for both V and B-V in this and 
in the following figures is 0.05mag. The histograms give the 
number of stars as a function of B-V: the total number is rep- 
resented by the full line histogram, while the dotted on is the 
Y=0.25 population and the dashed one is the population with 
variable Y 

not predict valus larger than Y t 0.36. We performed simula- 
tions with larger values for cr: while it is possible to achieve 
the end of the blue tail with Y<0.4, the mean RR Lyr period 
requires the presence of stars with Y up to ~ 0.37-0.38. These 
stars in fact have the (large) luminosity necessary to obtain the 
long observed periods. 

4.5. The number vs. helium distribution 

The number of stars as function of Y assumed in the simula- 
tions is given in the bottom panels of Fig. 6, 7 and 8. For our 
"best solution" we see that only 550/1452~38% of stars belong 
to the first stellar generation with "normal" helium Y=0.25. 
Anoher important feature is made evident by Fig.^| Stars with 
Y > 0.33 are necessary to achieve B-V < 0.5 mag, but a large 
fraction of these objects remains in the red clump, due to the 
large excursions in colour of helium rich structures (see Fig. 
O. Besides, the presence of structures with Y up to ~ 0.37 is 
required by the luminosity distribution of the clump itself (Fig. 
|6|l. This means that we must have a substantial amount of stars 
with Y > 0.33 in the cluster. In fact, in our "best solution" (§4.3 
and Fig. 5), we have -38% of stars with Y = 0.25, -33% of 
stars with 0.275< Y <0.33, and -29% of stars with Y > 0.33. 
About 15% of stars have Y > 0.35. In this last group are in- 
cluded the objects with Y = 0.39 - 0.40 (3.6%), which are the 
only ones to remain (almost) always at B-V < 0.0. 
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Fig. 6. For the simulation of Fig. 5, we show in the top panel the synthetic CM of the red clump region, with the details of 
the distribution in luminosity (line-dotted histogram, magenta). The continuous line gives the observed histogram, as in Fig. 
|2] Different symbols (and colours) denote different helium contents. Points: stars with Y=0.25 (red); open triangles (black): 
0.27<Y<0.32; asterisks (blue): 0.32<Y<0.38; stars (magenta): 0.38<Y<0.40. The dashed line indicates the approximate position 
of the RR Lyr variables, derived from Pritzl et al. 2003. The assumed age, the number of RR Lyr (RR) and the average period 
of the RR Lyr (P), and the fraction of HB stars bluer than B-V - 0.5 (R) in the simulation are listed at the top right. The lower 
panel shows the number distribution of HB stars vs. helium content corresponding to the CM diagram in the upper panel. The 
total number of the first stellar generation with Y=0.25 is given at the top right and labeled N(0.25). 



The exact numbers may be debated, as they depend in part 
on the age and mass loss along the RGB, and the extreme range 
Y = 0.39 - 0.40 can be avoided with a mass dispersion > 0.04 
Mq, but surely a substantial fraction (» 30%) has Y > 0.33. 

5. The R parameter and the bump location 

The ratio R of the number of HB stars to the number of red gi- 
ants above the HB level has been a powerful indicator of the he- 
lium content in cluster stars, as first suggested bv »Iben. (.1968a) . 
when the aim was to discriminate the helium mass fraction in 
the oldest stars of the Galaxy between values close to Y=0. 1 
or to Y=0.3. For NGC 6441, Salaris et al. (2004) find a value 
R=1.85 + 0.11, consistent with the theoretical expectation by 



the same authors for the cluster metallicity. On the contrary, 
the simulations bv'Sweigart & Catelan' d 19981) predicted values 
as large as 3.4 - 3.9, when they assumed main sequence helium 
abundances of Y=0.38 and Y=0.43 for all the cluster stars. For 
these values, the level of luminosity of the HB clump is highly 
enhanced and thus the number of red giants to be considered re- 
duced. This in fact is also predicted by the simple interpolation 
formula by Buzzoni et al. (1983) 

y = 0.380 log/? H- 0.176 (1) 

We do not explicitly compute the value of R from our sim- 
ulations, because it is easy to understand that we would find 
a "normal" value for R, due to the fact that the HB level in 
the R definition is the zero age value, which in this cluster is 
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Fig. 7. As in Fig.|6l but with no gap in the Y distribution (see 
the lower panel): it is not possible to avoid a peak in the lumi- 
nosity histogram substantially larger than the one observed 
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Fig. 8. As in Fig.|6l but with a gap of 0.035 between the Y of the 
red clump and the lower limit of the population with varying Y: 
two well defined peaks appear in the luminosity histogram. 



determined by the stars having the normally adopted value of 
Y=0.25. These are the lowest luminosity stars in the red clump 
and so the number of red giants to be considered remains the 
same as in the case with constant Y=0.25. In addition, there is 
only a small percentage of very hot stars with very high helium, 
whose HB lifetimes would exceed the Y=0.25 lifetimes, so that 
the correction to R due to t hese stars would turn out to be small. 
Therefore, the argument bv Lavden et al ] (11999') based on the R 
parameter can not be considered a proof against a high helium 
content as explanation of the RR Lyr long periods. 



iRaimondo et al notice that a uniform high helium 

abundance of the cluste r stars predicts a wrong locati on of the 
red giant branch bump (» Thomasi ll%7t llbeni fl968hl) with re- 
spect to the HB. Also in this case, our modeling is consistent 
with this latter parameter: a prominent RGB bump will occur at 
the luminosity predicted by models with Y=0.25, which is the 
abundance of ~40% of the stars. The stars having larger Y will 
produce smaller bumps, spread along the RGB at increasing lu- 
minosities, according to the number vs. Y distribution. Notice 
also that the prominence of the bump decreases with increas- 
ing Y. Further, the HB reference luminosity level in our simu- 
lations, as mentioned above, corresponds to the Y=0.25 level. 
Therefore, the quoted difficulty with a high helium population 
bv lRaimondo et alJ ( l2002l) does not apply to our model. 



6. Discussion 

The interpretation of the HB morphology of NGC 6441 in 
terms of a helium rich population has provided a coherent in- 
terpretation of the main features of the anomalous HB of this 
cluster, namely, the thick extension in luminosity of the red 
clump, the very long mean period of the RR Lyrs for a clus- 
ter of high metallicity, the extension into the blue of the HB. 
Modellization shows that a continuous distribution of helium 
content above Y=0.25 is not consistent with the luminosity dis- 
tribution of stars in the red clump, and that a helium disconti- 
nuity must be present, from Y=0.25 to a minimum value of 
Y^O.27, like in the cluster NGC 2808. We suggest that a more 
precise photometry of the HB would put into evidence a double 
peak in the luminosity distribution of the red clump stars. 

While it had been long suspected that the blue HB and 
the high RR Lyr luminosity in this cluster — and in its twin 
NGC 6388 — could be easily attributed to a high helium con- 
tent i n the small fraction of stars — 12%- constit uting the blue 
HB dSweigart & CatelanL Il998l ID' Antona & CaloL 2004), it 
comes as a surprise of these detailed simulations that a much 
more important fraction of stars has a very large Y: 29% have 
Y > 0.33 and ~14% have Y > 0.35. This latter figure is similar 
to the percentage of stars with Y ~ 0.4 found from the analy sis 
of the main sequence in NGC 2808 jP' Antona et al.'. 2005). In 
this respect, we conclude that both NGC 6441 and NGC 2808 
were able to form about the same fraction of very high helium 
stars, whatever their origin. 

The difference between the two clusters resides in the in- 
termediate helium population, which for NGC 2808 contains 
35% of stars, and covers o nly the range between 0.26 and 
0.29 dP' Antona et all l2005l) . In NGC 6441, this fraction is 
~ 48% and covers a heUum content from 0.27 to 0.35. Thus, 
any dynamical and chemical evolution model for these clus- 
ters must also explain the difference in this intermediate pop- 
ulation. We suggest that it is due to the diff'erent metallici- 
ties of the two clusters: in NGC 2808 (Z ~ 1 - 210"^) the 
maximum mass for AGB type evolution is smaller than in 
NGC 6441 (Z > 610"^). Therefore, the range of progenitor 
masses for the second stellar generation in NGC 6441 con- 
tains stars of larger mass which also have larger helium content 
(■Ventura. P' Antona. & MazzitemlEool) . 
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Fig. 9. As in Fig. |6l except that the standard deviation on V Fig. 10. The relation between B-V and heHum content Y for 
and B-V is reduced from 0.05 to 0.03 mag: the clump appears the case in Fig.|6l but a quite similar behaviour has been found 
separated in two groups. in all the simulations. 



One of the motivations to pursue the present work was to 
investigate whether an "extreme" high helium population is 
present in other clusters beside NGC 2808 and wCen. We found 
this extreme helium population in NGC 6441 (and this must 
also exist in the very similar cluster NGC 6388), and we also 
discovered that more than 60% of the cluster stars belong to the 
second stellar generation, having 0.27 t Y t 0.4. 

The interpretation of the HB morphology of NGC 6441 
makes it even more cogent to understand better which is the 
origin of these extreme helium rich stars. 
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